Diffusion induced decoherence of stored optical vortices 
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We study the coherence properties of optical vortices stored in atomic ensembles. In the presence 
of thermal diffusion, the topological nature of stored optical vortices is found not to guarantee 
slow decoherence. Instead the stored vortex state has decoherence surprisingly larger than the 
stored Gaussian mode. Generally, the less phase gradient, the more robust for stored coherence 
against diffusion. Furthermore, calculation of coherence factor shows that the center of stored vortex 
becomes completely incoherent once diffusion begins and, when reading laser is applied, the optical 
intensity at the center of the vortex becomes nonzero. Its implication for quantum information is 
discussed. Comparison of classical diffusion and quantum diffusion is also presented. 



Introduction Photons can carry orbital angular mo- 
mentum (OAM)jWhich can be created manipulated 
and detected @|. The OAM states are associated with 
vortices of a helical phase structure e %m ^ . Each vortex 
is a topological defect, characterized by a winding num- 
ber m, obtained from the 2nm phase twist around the 
vortex. The light with OAM is exemplified by the most 
commonly known Laguerre-Gaussian (LG) modes LG™ 
with p the number of nodes in radial direction [J, Q . En- 
tanglement through the OAM degree of freedom has been 
demonstrated for photon pairs from parametric down 
conversion Q. Since the winding number can be any 
integers, the OAM states have been proposed theoreti- 
cally for multiple-alphabet quantum cryptography with 
higher information density and a higher margin of secu- 
rity Other proposals, such as quantum coin tossing 
and violation of Bell inequalities, using the OAM states 
have been reviewed by Ref. @. 

One of the key challenges for optically based quan- 
tum information (QI) processing is the difficulty of stor- 
ing optical fields. It has been demonstrated [ljj, [0, G3 
that a superposition of the OAM states can be stored 
in a non-rotating BEC in terms of vortex states of the 
condensate. Meanwhile, the OAM states can also be 
stored in "hot" atomic ensembles using slow light tech- 
niques [H, G3, [U, GJ, [T3| - The information of pho- 
tonic states, namely the amplitude and phase, is con- 
tinuously transformed into Raman coherence, i.e., spin- 
density waves, of the atomic ensemble, which can be later 
retrieved. In practice, processes such as inhomogeneous 
magnetic field [l8| and/or thermal diffusion can lead to 
the decay of the Raman coherence. Then, to what extent 
can the optical vortex states be stored coherently? Are 
they going to be more robust than the Gaussian state? 
The answers to these questions will determine how the 
OAM arc used in QI. 

Generally speaking, the topological structure of vor- 
tices makes it a good candidate for QI [H, H(| Hl[ because 
it is stable against continuous deformations which cannot 
cause it to decay or to "unwind" . Actual studies of such 
robustness against various processes are of course neces- 
sary. In particular, one needs to study the robustness in 



the presence of diffusion [221 ] for the stored optical vortex 
state, which is crucial for applications such as quantum 
repeaters [23|, [24[ and multiple beam splitters for gener- 
ating entangled single photons 25|. Towards answering 
this question, Pugatch, et al. [l7| showed that in the 
presence of diffusion, the dark center of a stored OAM 
mode is well preserved, and the dark center of a Gaussian 
mode generated by blocking its center disappears after a 
short time of diffusion. They attributed the stable dark 
center of the stored OAM mode to the robustness of the 
topological nature of vortex in the presence diffusion. 

In this Letter, we provide a careful study for the ro- 
bustness of the stored vortex states (exemplified by p = 
unless otherwise stated) in the presence of diffusion. We 
find that the stable dark center of the vortex states is 
not directly associated with the topological robustness. 
And the vortex states are actually more vulnerable to 
diffusion than the Gaussian state. This is because: (1) 
the diffusion is a global process during the storage and 
it can destroy the topological order; (2) the readout pro- 
cess only measures a specific output modes corresponding 
to Raman coherence but not the full state; (3) for vor- 
tex states, the Raman coherence interferes destructively, 
while the full quantum state evolves differently. 

Our basic assumption is that there is no spin exchange 
or interaction between atoms, the evolution of coherence 
obeys diffusion equation p = D\7 2 p, where p and the 
D are the density matrix and the diffusion coefficient, 
respectively. We consider a three-level lambda system as 
generally used for light storage [H, [3| : the weak probe 
laser is applied between the ground state |1) and the 
excited state |3), and the pump laser is addressing the 
transition between |3) and another ground state |2). 

Comparison between the decay for vortex states and 
vortex-free states Using the propagator of the diffusion 
equation, the expected atomic coherence of a LG mode 
after diffusion of duration t is [l7| 

p u (r,t) = { ~ 9/n) A m (r, ^tjw )e~ tme (I) 



with g the vacuum Rabi frequency for the probe field, 
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il the Rabi frequency for the pump laser, s(t) = 
(wq + 4Dt)/uiQ an evolution factor, wo the waist, 
and the radial profile of a LG mode A m (r, wq) = 
{l/w t) ) v /2P/Trml(V2r/w ) m exp(-r 2 /w$) where P is 
the total intensity. Note that Eq. fl} applies to the Gaus- 
sian mode (m = 0) as well. For example, the evolution 
of stored coherence is shown in Fig. []Ji and Fig. [5Ji for 
m = 1 and m = 0, respectively. At large rs, the coher- 
ence is homogeneously zero due to the exponential factor 
in A m . Note that Eq. JT]) describes both a spread of the 
coherence, indicated by \J s(t) inside A m , and a decay of 
the coherence compared a purely coherent spread. The 
decay is given by \J s(t) m+1 in the denominator. Then 
F = l/s(t) m+1 gives the fidelity of the stored coherence. 
F shows that different angular momentum states ms have 
different coherence decay factors. Although we will come 
back to the decay of coherence later, we note here that the 
larger OAM m, the larger the decay factor. The Gaus- 
sian mode, which has no phase singularity, has smaller 
coherence decay factor than all vortex states! This is be- 
cause without vortex, the coherence is always in phase 
for different locations and there is no destructive inter- 
ference to destroy the coherence. One can show that the 
retrieval efficiency along forward direction is the same as 
the fidelity F. Therefore, we come to a counter-intuitive 
result — a Gaussian state has a higher storage fidelity 
(or retrieval efficiency) than the LG modes. 

This may make the applications of optical vortex (with 
storage) for QI questionable. Even though the excita- 
tion loss of stored Gaussian or optical vortex state can 
be identified as a detected error and the threshold for 
detected errors allowed by quantum computation are re- 
alistic [2H | , fast quantum computation with reasonable 
resource overhead still requires relatively low error rates. 
This is also true for quantum communication [23l , [24j], 
where, although the excitation loss can be fully controlled 
by its intrinsic purification, a fast quantum repeater still 
needs a fairly high retrieval efficiency. 

The scaling factor y/ s(t) in A rn of Eq. |T]) means that 
the functional forms of coherence are preserved for both 
the LG mode and the Gaussian mode. This implies the 
functional form stability of stored coherence against dif- 
fusion does not require topological defects. Indeed, the 
disappearance of the dark center of the center-blocked 
Gaussian mode after diffusion [l7| demonstrates as the 
stability of the Gaussian mode: diffusion tries to restore 
the nonzero intensity at its center. Of course, the restora- 
tion can also be understood by decomposing the center- 
blocked Gaussian mode to LG° and noting that p ^ 
modes decay faster than p = mode and what is left 
after some time of diffusion is just the Gaussian mode. 

From the above discussion, it is clear that for the LG 
light (m 7^ 0), the coherence at the center r — stays zero 
during diffusion. However, as we show now, population 
at r = does not stay zero. To simplify the discussion, 
we assume, at time t = 0, the populations of atoms are 
Pll (r,t = 0) = 1 and p 22 {r,t = 0) = \ Pl2 (r,t = 0)| 2 , 
which is a good assumption for strong pump and weak 



probe lasers as usually used in light-storage experiments. 
At time t after diffusion, we have pii(r, t) = 1 and 



4e * Dt+ < P (32DH 2 + r 2 wl + 4Dtw$) fS 
#22 M) = jzz: ^3 ( 2 ) 

for m = 1 vortex state and 



2e 8Dt+ ™5P 
ir{8Dt + w 2 ) 



(3) 



for a Gaussian mode. The evolutions of populations are 
plotted in Fig. [lfo and Fig. [2}x While it seems that co- 
herence only diffuses outwards in Fig. [T^,, Eq. @ and 
Fig. [lb clearly show that diffusion goes in all directions 
as it should be. The outwards moving coherence during 
diffusion is because the interference cancels the inwards 
diffusing coherence. In contrast, the population does not 
interfere with each other and thus the diffusion towards 
the center is clearly seen. Indeed, the population at the 
center quickly approaches a global maximum as time in- 
creases (Fig. [TJd). We also note that integration of p 22 
over the whole space is conserved during the diffusion. 

Diffusion induced decoherence We have seen that diffu- 
sion for coherence and population obeys different equa- 
tions (Eq. |T][3])). This brings phase decoherence for the 
stored coherence. To characterize the decoherence, we 
define a coherence factor f = Pi "^ + . v , where 77 is an in- 

. PbbPcc+V 1 _ ' 

finitesimal number. / is a linear function with / = 1 for a 
pure state and / = for a completely mixed state. Thus, 
/ is a good parameter to describe the (local) coherent 
property. As a specific example, we plot the coherence 
factor / in Fig. QJ: and Fig. [2t. Figure [IJ shows that 
right after diffusion begins, coherence factor f(r = 0) of 
the stored vortex drops from 1 to because at r = 0, 
the population becomes nonzero when diffusion starts 
while its coherence keeps zero. We note that such sud- 
den changes within a infinitesimal time are very uncom- 
mon in physical processes. As diffusion time increases, 
f(r) drops faster for rs that do not have much initial 
inhomogeneous coherence, and f(r) approaches zero at 
very large time. This latter result holds for a Gaussian 
mode as well (shown in Fig. [2}:). We also note that at 
rs where the diffusion of stored inhomogenous coherence 
and population has not arrived at, the coherence factor / 
stays at 1 because all population is in |1) and it is a pure 
state. But the weight p 22 , justified by its non- negative 
and conserved integration, of these coherence factors for 
the retrieved light approaches zero. 

Here is how coherence factor / may be obtained from 
the experiments. When the reading pulse is applied, the 
incoherent part will be retrieved as fluorescence in all di- 
rections. Collecting both the fluorescence and the coher- 
ent emission then allows extracting /. Of course, setting 
the detector at the forward direction as generally used, 
e.g., in [T3| , can only collect a finite ratio of fluorescence, 
while the coherent part is collected by the detector [27j . 
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FIG. 1: (Color online) Effects of diffusion on a stored vortex 
with topological charge m = 1. Plotted are coherence pi2 (a), 
population P22 (b) and coherence factor / (c) as functions of 
radius, for different diffusion times, wo, D are the waist of 
the LG mode and the diffusion coefficient, respectively. 



FIG. 2: (Color online) Effects of diffusion on a stored Gaus- 
sian mode. Plotted are coherence pn (a), population P22 (b) 
and coherence factor / (c) as functions of radius, for different 
diffusion times, wo, D are the waist of the Gaussian mode 
and the diffusion coefficient, respectively. 



What we want to emphasize is that incoherence makes 
intensity at the center of a retrieved vortex nonzero, pro- 
hibited by a coherent vortex state. As diffusion time in- 
creases, nonzero intensity at the center builds up. There- 
fore, generation of mixed state makes additional loss of 
retrieved fidelity compared with vortex-free case. Diffu- 
sion of the population makes visibility decrease and fi- 
nally kill the vortex. Of course, this part of reduction of 
fidelity can be alleviated by using spatial filtering of the 
optical mode to prevent spontaneously emitted photons 
from going into the detectors. In this case, the retrieval 
probability is just the fidelity F = l/s{t) m+1 . 

The different collection efficiency of fluorescence pho- 
tons and coherent photons by a forward set detector helps 
explain why the visibility of the center-blocked Gaussian 
mode disappears very quickly, while the hole of vortex 
disappears very slowly [17] . The homogeneous phase of 
the center-blocked Gaussian mode makes nonzero coher- 
ence inside the hole after diffusion and thus the disap- 



pearance of hole once the coherence is read. This is in 
contrast with a stored vortex, whose coherence at the 
center remains zero all the time. Its nonzero intensity 
at the center comes only from incoherence of the center. 
However, were most of fluorescence photons collected by 
the detector, the dark center of a vortex would disappear 
quicker than that of the Gaussian mode. It is the spatial 
filtering of the optical mode that helps to overcome the 
fluorescence from the center. 

Dependence of decay on order of phase singularity We 
now come back to the decay of coherence. We have noted 
that the coherence of stored LG modes decays according 
to a power law F = l/s(t) m+1 . The larger the order of 
phase singularity, given by to, is, the faster the coherence 
decays. An additional example of exponential decay rate 
2Dk 2 due to diffusion of a plane wave e~' lkx , which is 
faster than the power law decay, also corroborates this 
idea, because a plane wave has infinity number of phase 
singularities. Incidcntly, the limit of F for large to does 
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not go to exponential decay of a plane wave is because the 
amplitude of coherence in LG modes is not a constant, 
different from the case of a plane wave. We also note that 
the larger k of a plane wave, the larger decay rate, which 
is because larger k means the pattern have higher spatial 
frequency and the diffusion becomes easier to cancel the 
coherence. Such diffusion of a plane wave happens if the 
pump and probe lasers have different wave vectors. 

A few more remarks on the different decay behaviors 
arc in order. First, as a direct application of our discus- 
sion, the less phase gradient, the better for stored coher- 
ence against diffusion. For example, diffusion of stored 
general LG™ with both winding number and number of 
radial nodes being non-zero @ induces faster decay for 
p > than p — 0. Furthermore, we note that although 
the number of nodes in the coherence does not change 
with diffusion time, the positions for the off-center ra- 
dial nodes change, which is different from the center one. 
The non-moving position of the center node comes from 
geometric symmetry of the vortex. Since the decay rate 
depends on both |m| and p for the mode LG™, if OAM 
states are to be used as bases for quantum information, a 
preferred basis to reduce loss of entanglement for quan- 
tum information is actually two modes with same p but 
opposite m. Second, since diffusion induced decay of 
the total intensity is not exponential for Gaussian and 
LG modes, the exponential decay at a rate 20, OOOs^" 1 of 
Ref. [TtJ is not determined by diffusion. 

The classical diffusion vs the quantum diffusion So far, 
we have only discussed decoherence induced by the clas- 
sical diffusion associated with the inhomogeneous dis- 
tribution of coherence. We note that decoherence can 
also happen in a homogeneous system as a result of the 
quantum diffusion. The quantum diffusion in light stor- 
age system happens when pump and probe lasers cou- 
ple different momentum states, which introduces deco- 
herence (28|. But this decoherence is reversible, say by 



photon echo techniques, in contrast with the classical 
diffusion. This is because the quantum diffusion is de- 
scribed by the complex Schrodinger equation with i in it, 
while classical diffusion does not. Finally, we note that 
when the same momentum states are coupled by choosing 
pump and probe lasers to have the same wave vectors, 
the quantum diffusion disappears [28| . fncidently, inho- 
mogenous magnetic fields also induce quantum diffusion. 

Discussion Our results indicate that diffusion actually 
introduces more decoherence in a stored vortex mode 
than a stored Gaussian mode, which may have impor- 
tant implication for quantum information. Of course we, 
however, do not rule out that in other processes, such as a 
quantum gate operation, vortex states are possibly much 
better than the topological-free state. Nor did we dis- 
cuss diffusion-free systems such as BEC [29[ and bound 
excitons in semiconductors [30l . l3ll ] . 

Conclusion We found that during diffusion, the co- 
herence of stored vortex states decays faster than that 
of Gaussian states. This is surprising because vortex 
states are associated with topological properties, and are 
presumably considered as more stable than topological- 
defect-free states. The underlying reason is that diffusion 
is a non-local process. More generally, the less phase 
gradient of a stored coherence, the better for it against 
diffusion. Furthermore, calculation of coherence factor 
showed that the center of stored vortex becomes com- 
pletely incoherent once diffusion begins, and when read- 
ing laser is applied, the optical intensity at the center of 
the vortex builds up. It's implication for quantum infor- 
mation was discussed. Finally, we compared the classical 
diffusion and the quantum diffusion. 
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